Coherence of resonant-tunneling transport in terahertz quantum-cascade lasers 
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We develop simple density-matrix models to describe the role of coherence in resonant-tunneling 
(RT) transport of quantum-cascade lasers (QCLs). Specifically, we investigate the effects of coherent 
coupling between the lasing levels with other levels on the transport properties and gain spectra. 
In the first part of the paper, we use a three-level density-matrix model to obtain useful analytical 
expressions for current transport through the injector barrier in a QCL. An expression for the slope 
discontinuity in the current-voltage characteristics at the lasing threshold is derived. This value is 
shown to be a direct measure of the population inversion at threshold, and contradicts the previously 
held belief of it being indicative of ratio of the laser level lifetimes. In the second part of the paper, 
we use density matrices to compute the gain spectrum for a resonant-phonon terahertz QCL design. 
The large anticrossing of the doublet of lower radiative levels is reflected in a broad gain linewidth 
due to a coherent RT assisted depopulation process. At certain bias conditions, the gain spectrum 
exhibits double peaks which is supported by experimental observations. 
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I. INTRODUCTION 

The role of sequential resonant-tunneling (RT) in semi- 
conductor quantum-wells for charge transport [1] as well 
as light amplification [2J has been discussed since the 
early 1970s. Subsequently, multiple proposals were in- 
troduced in the 1980s to predict intersubband lasing 
in semiconductor superlattice structures. With the ad- 
vent of a unique RT based charge injection scheme, a 
unipolar mid-infrared (IR) intersubband laser, the so- 
called quantum-cascade laser (QCL), was first demon- 
strated in 1994 Mid-IR QCLs have since under- 
gone rapid progress and room-temperature continuous- 
wave operation is now obtained over the A ~ 4 — 10 /im 
spectrum. Based on similar design principles, the op- 
crating frequency range of these devices has also been 
extended to terahertz frequencies (y ~ 1 — 5 THz, 
A - 60 - 300 /zm) Q. However, due to additional chal- 
lenges associated with the charge transport in these low- 
frequency devices, terahertz QCLs are still required to be 
cryogenically cooled [E B 0] • ^ then becomes increas- 
ingly important to better understand the transport pro- 
cesses in terahertz QCLs to improve the existing designs, 
with a goal to ultimately achieve room-temperature op- 
eration. 

The best temperature performance for terahertz QCLs 
is obtained with the resonant-phonon (RP) designs M,M , 
and operation up to 186 K without magnetic field [lfj, 
and 225 K with magnetic field [ll[ has been demon- 
strated. For QCLs, almost universally across all different 
designs and operating frequencies, RT plays an important 
role in the electrical injection process to populate the up- 
per radiative level [12, [HI- For the resonant-phonon tera- 
hertz QCLs (RPTQCLs), RT also plays an important role 
in the lower level depopulation. Newer designs with RT 
based extraction have also been demonstrated recently 
achieving very low threshold current densities 0, fl5f . 
The radiative levels in a terahertz QCL have small en- 



ergy separation in the range of hw ~ 4 — 20 meV 
{y ~ 1 — 5 THz). The potential barriers in the mul- 
tiple quantum-well (QW) structure are therefore kept 
thick to obtain tightly-coupled levels in order to maintain 
selectivity of the injection as well as the depopulation 
processes. This makes the coupling energies (as charac- 
terized by the anticrossing energy 2hVLij for the coupled 
levels i and j) similar in value to the low-temperature 
energy level broadening, which is expected to be of the 
order of a few meV. Consequently, the loss of coherence 
(or dephasing) in the RT process has a significant bearing 
on the electron transport across the barriers [ll|. 



In this paper, simple density-matrix models are used 
to incorporate the important role of coherence in the RT 
transport of QCLs. We consider the case of tight binding 
for intermodule transport, since the intermodule energy 
anticrossings 2Mlij are typically much smaller than the 
intramodule energy level separations, where each mod- 
ule can consist of one or more QWs and the QCL struc- 
ture is a periodic repetition of one or more such mod- 
ules. In the first part of the paper, we use a density- 
matrix model similar to the 2-level model first proposed 
by Kazarinov and Suris [r7| and then later also discussed 
in Refs. [H, EH [H among others. We modify the model 
to extend to 3 levels to describe a QCL in a more general 
sense, and derive analytical expressions for current trans- 
port across the injector barrier both below and above the 
lasing threshold. Despite the simplicity of the model, we 
can gain useful information about QCL operation from 
the derived results that includes information about the 
possible temperature degradation mechanisms in RPTQ- 
CLs. The results presented in this section are applicable 
to mid-IR QCLs as well. In the second part of the paper, 
we use a similar model to estimate the optical gain spec- 
trum of RPTQCLs. A coherent RT assisted depopulation 
process is stipulated to be the cause of broad gain band- 
widths that are typically observed in these lasers [2p] |. 
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FIG. 1: Plot of the magnitude squared envelope wavefunc- 
tions for a 3-level QCL design with two different sets of basis 
functions. Plot (a) is for the "extended" scheme, where energy 
splitting due to the injector anticrossing 2hQ 1 i 3 is visible, and 
the wavefunctions are calculated for a potential profile // ex t 
as it appears in the figure. Plot (b) is for the "tightbind- 
ing" scheme, where the potential i/ t b is formed by making 
the barriers at the boundaries of a module sufficiently thick 
to confine the wavefunctions within the module. 



II. CURRENT TRANSPORT IN A 3-LEVEL QCL 

An accurate estimation of the electrical transport char- 
acteristics of a QCL is quite challenging, and more so 
for a terahertz QCL. Computationally intensive numeri- 
cal techniques based on Monte-Carlo simulations [l6j] and 
nonequilibrium Green's functions [21] have produced rea- 
sonably good results; however, an accurate description 
of the temperature degradation mechanisms in terahertz 
QCLs is still lacking. Also, such computational tech- 
niques fall short of providing a good intuitive picture of 
QCL operation. Hence, more often than not, analysis 
based on simple rate equations are used during the de- 
sign process. The 3-level model presented in this section 
is a step toward developing a simpler transport model 
that still captures the nuances of QCL operation and 
also provides a better understanding of how various pa- 
rameters affect its characteristics both below and above 
the lasing threshold. 

For simplicity yet still captures the essentials, in what 
follows we discuss a 3-level QCL design with two wells per 
module as shown in Fig. Q] [2^, , in which level V is the 
injector state, and the radiative transition is from level 
3^2. Even though typically QCL designs include many 
more levels per period, the results obtained here hold for 
the triplet of levels consisting of the injector level and the 
two radiative levels (levels 1', 3 and 2 in Fig. [TJd), regard- 
less of a particular design, and including the mid-IR de- 
signs. Due to the various intrasubband and intersubband 



scattering mechanisms that contribute to quantum trans- 
port, the electron wavepackets undergo dephasing, which 
manifests itself as broadening of the energy levels. The 
tunnel-coupling between the levels of the same module 
is strong, therefore, dephasing has a relatively negligible 
effect on the intramodule intersubband transport fig . 
However, the intermodule transport across the injector 
barrier is critically affected by the dephasing time of the 
wavepackets since the 1' — > 3 tunneling time (half of the 
inter-well Rabi oscillation period Tr/Qva) is now longer 
and of the same order. To incorporate the role of coher- 
ence in the 1' — > 3 RT transport, we use density-matrix 
(DM) rate equations corresponding to the tightbinding 
basis states of Fig. QJb) [lq . The time evolution of the 
3x3 DM for this basis set can be written as 12211 



dt 



P(l',3,2) 



H, 
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where, H includes various radiative and non-radiative 
perturbation potentials to the conduction-band potential 
H cx t- The advantage of writing the DMs with the tight- 
binding basis of Fig.[ljb) is two-fold. First, it allows in- 
clusion of dephasing phenomenologically by means of de- 
cay of the coherences associated with levels V and 3 due 
to non-zero off-diagonal terms in H ext corresponding to 
those levels. Second, whereas the eigenstates of i? cxt de- 
pend sensitively on the externally applied electrical bias 
because of coupling between spatially separated states, 
the chosen basis states remain more or less independent 
of the bias in the range of interest (i. e. close to the "de- 
sign bias" corresponding to the 1' — 3 alignment), which 
keeps the form-factors for various intramodule scattering 
rate calculations approximately invariant of bias. In a 
first order approximation, we assume RT to be indepen- 
dent of the electron wavevector in the plane of the QWs. 
The subband lifetimes are assumed to be averaged over 
the electron distribution within a subband and are to be 
calculated within the Fermi's golden rule approximation. 
Expanding equation fT]). we then obtain 
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where pu (= n{) is taken as the number of electrons per 
module in level i and pij is the coherence (also known 
as polarization) term for levels i, j. In H ex t, the diagonal 
terms are the level energies given by H nn = (n\H cxt \n) w 
E n and the level anticrossings are represented in the off- 
diagonal terms as H mn = (m\H eKt - H th \n) w —Ml mn . 
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For the lifetimes, r^- is the intersubband scattering time 
from i —y j, -J— = — h ?r — h ^ is the dephasing rate 

for the coherence term pij that consists of lifetime broad- 
ening terms as well as a phenomenological broadening 
ter m Ty due to interface-roughness and impurity scatter- 
ing [l6|, and r s t is due to the radiative stimulated emis- 
sion above the lasing threshold (— > oo below threshold). 
The backscattering times t 2 3 and r 12 can also be included 
in equation @ if relevant for a particular design. Note 
that f^32 = for the chosen bases since subbands 3 and 
2 are the eigenstates of the tightbinding potential H t b hi 
Fig.QJb). By such a choice, the role of coherence in 3 — > 2 
transport cannot be included. This does not introduce a 
large error in estimating the current flow since 3 and 2 are 
strongly coupled. In other words, if tightbinding bases 
were chosen for the radiative subbands 3 and 2 across 
the middle "radiative" barrier of the two-well module, 
f2 32 would be large and the Rabi oscillation period ir/fl^ 
would be much smaller than the dephasing time T|| 3 2 for 
the 3 — > 2 tunneling, making t"|| 32 inconsequential for the 
3^2 transport. 

As a simpler alternative to more advanced methods, 
equation |2j) can model some aspects of QCL transport 
fairly accurately. Although it could be solved numeri- 
cally, we seek analytical expressions that could provide 
a greater understanding about the effect of various pa- 
rameters on the transport. Toward that goal, we assume 
VL\'2 ~ akin to a unity 1' — > 3 injection selectivity. This 
assumption is in general valid for mid-IR designs due to 
a large radiative level separation E32 (= £3 — E 2 ), and 



also becomes reasonable for the diagonal terahertz de- 
signs [lCj |. Within this approximation, equation @ can 
be reduced to that with 2x2 matrices as follows 
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P33-P22 f rom equation ((2]) is substituted by 
— in equation (J3J, which holds when Oj/2 = as 
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a statement of current continuity. Additionally, we can 
write the following equation for p 2 2 below and above the 
lasing threshold 
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where, Anth = (P33 — /022)th is the population inversion 
at threshold that is assumed to remain constant beyond 
threshold. Equations © and (g]) can be solved ana- 
lytically for steady-state (4t — » 0). With the constraint 
"tot = (pn + P22 + P33), where n to t is the total number 
of electrons per module and is a constant, the following 
expressions are obtained for the current / [22j | 
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where, e is the unit charge and hAy 3 = Ey — E3 is 
the energy detuning between levels 1' and 3 that is a 
function of the externally applied electrical bias. For the 
case of ultra-short lifetime of the lower radiative state 
(t"21 *C 2T32), the below-threshold term in equation (O 
becomes 

j = [ejntot / 4 »i'3 T ll r 3 \ 

2r 3 y^l^Ta + Al^ + l) 

. . . (T21<2T32, -f<-fth) (6) 

which reproduces the commonly used expression in liter- 
ature UllI. 

Since T21 <C 2T32 typically holds for a QCL design, 



the below-threshold term for the 3-level model in equa- 
tion (O does not provide any new insight from equa- 
tion ©• However, the advantage of the model lies in 
describing transport above threshold. It is instructive 
to write expressions for maximum current (I max ), which 
flows at the 1' — 3 resonance (Ai'3 = 0) for the case 
of a coherent RT process, i. e. a large injector coupling 
[IllHil]- The following expressions are obtained for 
I max and the level populations when threshold could not 
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be attained 
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The results in equation (O are slightly different than 
that obtained using a two-level model in Ref. [H as a 
consequence of adding an additional level (i. e. the lower 
laser subband 2) in these calculations. However, the com- 
monly used two- level model result 7 max = |e|n to t/(2r3) is 
recovered in the limit tix -C T32 as it is to be expected. 
Equation (|Tj) suggests that J max for coherent injection is 
independent of fii'3, and therefore the thickness of the 
injector barrier, and is limited by the upper state lifetime 
T3. Typical values of £lv 3 that bring a QCL within this 
limit could be determined from Fig. [2ja). The value of 
Imaxj however, is different if enough population inversion 
could be attained prior to the 1' — 3 resonance and the 
device starts lasing, in which case stimulated emission 
lowers the upper state lifetime. In such a scenario, the 
maximum current is limited by the lower state lifetime 
T21 instead, and is given by 
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Equation (JSJ) implies that a QCL with a large injector 
coupling and a short lower state lifetime T21 can 
possibly obtain a large dynamic range in current due to 
a large I max (see Figs. [3]; and|3jl). However, it may be 
noted that for coherent injection, the gain linewidth of 
the QCL becomes additionally broadened as the injector 
coupling r2i'3 is increased. Consequently, a larger Anth 
is required to attain a certain value of peak gain, which 
mitigates the increase in the value of I max . This will 
become more clear in section Infl from the gain spectrum 
calculations for some specific cases. 

Approaching threshold, the value of current that must 
flow through the QCL structure to establish a population 
inversion of Anth is calculated to be 
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While the above expression for 7 t h can also be derived 
from a simple rate equation analysis, the minimum value 



of O1/3 needed for a current /th to flow through the struc- 
ture is derived from the present 3-level DM model as 
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This result is obtained from equation as a necessary 
condition for the two expressions to have the same value 
for some particular bias A]^. An estimate of the min- 
imum value of f2]/3 required to meet this condition can 
be obtained from the plots in Fig. [2](b). Figure [2jb) can 
also be interpreted in a different way, since it determines 
the maximum value of population inversion that can be 
attained for a given Qys. For example, in the limit 
T21 — > 0, a maximum of An = n to t/2 can be obtained 
if 40?, 3 t||T3 > 1. 



A. Current-voltage characteristics for typical 
parameters 

To gauge the effect of various parameters on QCL's 
performance, Fig. [3] shows calculation of the current- 
voltage (I-V) and differential-resistance-voltage (TZ-V, 
1Z = 4j) characteristics for a range of parameters us- 
ing the analytical expressions in equation (JSJ. The I-Vs 
show a discontinuity in slope at the occurrence of the las- 
ing threshold, as was hrst discussed in detail in Ref. H2I 
However, any discussion about the differential resistance 
in the present context is left until later. For the I-Vs the 
main features to be noted are the values of the threshold 
current density J t h and the maximum current density 
Jmax- In general for any design, the goal is to obtain 
a larger dynamic range for lasing J max /Jth while also 
keeping the value of J t h low. A larger dynamic range 
leads to greater amount of the optical power output, and 
also a higher operating temperature T max . As it is to 
be expected, the plots show that the laser performance 
improves by making T31 and T32 larger, and r 2 i smaller. 

We now compare the 3-level model I-Vs in Fig. [3] 
to those experimentally measured for two different RP- 
TQCL designs. Fig. [5] shows typical results from a 4- 
level [13, Hi H3 and a 5-level @, Q RPTQCL respec- 
tively. The structures in Fig.[D(a) andQJb) can be qual- 
itatively analyzed with the 3-level model of this section 
by using the following expression for the effective lifetime 
of the lower laser level 2, which is depopulated through 
RT into 2a, and subsequently via electron-longitudinal- 
optical (e-LO) phonon scattering into the injector level(s) 
(i. e. the resonant-phonon scheme). The following ex- 
pression can be derived using a 2-level density-matrix for 
levels 2 and 2a exclusively and assuming 2 — > 2a tunnel- 
ing to be the only mechanism for carrier injection into 
level 2a [H 
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FIG. 2: (a) A plot of the factor Aflp^Ts to determine the regime of RT transport through the injector barrier (^> 1 coherent, 
<C 1 incoherent). A phenomenological value of T 2 * = 0.75 ps is assumed, (b) A related plot to determine the maximum value of 
population inversion (713 — 712) that can be attained for a given fii'3, calculated using equation (|10[1 . The top axis converts An 
on the bottom axis to corresponding values of peak gain for the following typical parameters: ntot /Volume = 5 x 10 lj cm -3 , 

fs2 = 0.5, Av = 1.0 THz, and using the expression g ~ 70 ^"^^/(^thz)? '-^ [cm -1 ] for GaAs material [f| that assumes a 
Lorentzian gain linewidth of Av (note that f$2 is the oscillator strength of the radiative transition). 
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FIG. 3: Calculations of the current density (//Area) as a function of the energy detuning (which is proportional 

to the external bias voltage V, as in equation I12[). and its inverse slope j§y— ^j- 3 - (which is proportional to the differential 
resistance 1Z, as in equation I13fl . The calculations are done for a range of parameters using equation ([SJ for the 3-level model. 
Following default values are chosen corresponding to the typical values for RPTQCL designs: 2fi,fi 1 / 3 = 1.5 meV, T 2 * = 0.75 ps, 
ffth = 40/cm [2^|, T21 = 0.5 ps, T32 = 3.0 ps, T31 = 5.0 ps. A doping density n to t /Volume = 5 x 10 15 cm~ 3 is used. The 
radiative transition is assumed to have a Lorentzian linewidth of Av = 1.5 THz (the broad linewidth is based on the findings of 
section iHlll and an oscillator strength of jz2 = 0.6 to determine the population inversion Ann, required to attain a particular 
value of threshold gain g^. The thin (red) portion of the curves corresponds to I < J t h, while the thick (blue) region is for 
I > Ith where the occurrence of threshold is marked by a discontinuity in 1Z. Each of the subplots show the variation of the 
I-Vs and the IZ-Vs when only a single parameter is changed, the others being kept at the values mentioned above. Note that 
additional linewidth broadening due to coherent effects (discussed in section |HT)| is not considered for calculations of the curves 
above threshold. 
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(b) 5-level resonant-phonon design 
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FIG. 4: Typical experimental I-Vs versus temperature for 
(a) a 4-level, and (b) a 5-level resonant-phonon terahertz 
QCL, respectively. The corresponding band diagrams show 
tightbinding wavefunctions at design-bias calculated by split- 
ting the QCL module into multiple submodules across the 
relevant barriers. The radiative transition is from 3^2, 
and the depopulation of the lower level is via 2 — » 2a RT 
and 2a — > 1/la electron-longitudinal-optical phonon scatter- 
ing where £21 ~ hu^o. The I-Vs are measured in continuous- 
wave (cw) operation for metal-metal ridge lasers. The plot in 
(a) is from a 3.9 THz QCL labeled OWI222G [3 of dimen- 
sions 60 )im x 310 fim (T maXiC w ~ 76 K), and that in (b) is 
from a 2.7 THz QCL labeled' FL178C-M10 of dimensions 
35 |iim x 670 /im (T max ,cw ~ 108 K). The onset of lasing re- 
sults in a change in the slope of the I-V, which is indicated 
by arrows for the curves recorded below the T max , cw . 



The effective lower-level lifetime T2 iC s, which is also 
equivalent to the tunneling time through the collector 
barrier, is plotted for three different values of the collec- 
tor anticrossing 1MXi,2a in Fig- El The best performing 
RPTQCLs have large collector anticrossings in the range 
of 2hD,2.2a ~ 4 — 5 meV, in which case T2 yC s varies little 
in the operating bias range close to the 2 — 2a resonance 
(h/S.2 : 2a = 0). For example, 2hVl2 : 2a ~ 4.7 meV for the 
4-level QCL in Ref. [l0| and the energy detuning SA 2j 2a 
approximately varies from —2 meV to 2.5 meV in oper- 
ating bias range of the design. As seen from Fig. EJa), 
T2,eff changes little within such a range of detuning. 

We can similarly take the combined population n to t = 
ny +712 + ^3 to be approximately bias independent. The 
population of level 2a is likely to be negligible due to its 
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FIG. 5: Plots of the effective lower-level lifetime T2, c s (also 
equivalent to the tunneling time through the collector barrier) 
for the resonant-phonon designs of the type shown in Fig. [4] 
calculated using equation (|lip . Plot (a) shows T2. c fr versus 
the energy detuning ftA2,2o between the levels 2 and 2a and 
(b) shows T2, c ft versus the phenomenological dephasing time 
T2 calculated at the 2 — 2a resonance bias (/iA2,2a = 0). 
T2a ~ 0.2 ps is assumed owing to the typical value obtained 
in GaAs/Alo.isGao.ssAs quantum wells for E% a ,\ ~ ft^LO- 



very short lifetime T2 a ~ 0.2 ps, and that of level la in 
Fig. IHb) is likely to vary negligibly with bias in the oper- 
ating bias range of interest due to the typically large val- 
ues of the intra-injector anticrossing 2Ml\ a i ^4 — 5 meV 
in the 5-level design Fig. 0Jb) . Given the bias indepen- 
dence of ntot = Tiy + ri2 + 713 and the lower-level life- 
time T2 : eff for the resonant-phonon designs in Fig. [4] we 
can now qualitatively compare the 3-level model I-Vs in 
Fig. [3] to the experimental continuous-wave (cw) I-Vs 
in Fig. SJ For the latter, it can be noticed that J max , 
which is characterized by the occurrence of a negative- 
differential resistance (NDR) region in the I-V, changes 
little up until the maximum operating temperature of 
the laser T roaX)Cw . However, beyond r max ,cw, the value of 
J mas increases with temperature. This characteristic is 
best described by the calculated I-Vs in Fig. Eta) since 
the above-threshold current expression in equation (JS]) 
is independent of the upper state to lower state lifetime 
T32 to a first order. Hence, it is likely that one of the 
dominant temperature degradation mechanism for RP- 
TQCLs is a reduction in the upper state (it) to lower 
state (I) lifetime, which is mostly attributed to the ther- 
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mally activated LO phonon scattering from the upper 
state [5, 6]. The recently developed terahertz QCL de- 
sign with a very diagonal radiative transition is based 
on this observation to effectively increase the u — > I life- 
time at high temperatures [Io( |. Additionally, it seems 
less likely that an increase in the effective lifetime of the 
lower state (potentially due to 1 — > 2 thermal backfill- 
ing or increased dephasing for the 2 — > 2a RT that af- 
fects T2 : eff as in Fig. [5p), or a broadening of the gain 
linewidth [28j are more dominant temperature degrada- 
tion mechanisms than that mentioned above; otherwise, 
the observed experimental behavior would have corre- 
sponded more closely with Figs.[3^e) or[3lc), respectively. 
That is, the peak current density J ma x would decrease at 
elevated temperatures. 




Bias (V) 

FIG. 6: (a) Conduction band diagram of a recently demon- 
strated two- well terahertz QCL design, labeled TW246 [23| . 
Only three levels (1,2 and 3) participate in transport at 
low temperatures; however, at higher temperatures additional 
parasitic levels (4 and 5) are also likely to contribute to the 
current flow, (b) Experimental cw I-Vs and IZ-Vs at different 
heat-sink temperatures from a 20 /im x 1.56 mm metal-metal 
ridge laser with a T max ,cw ~ 63 K. The lasing threshold is 
marked by a discontinuity in the TZ curves that is shown more 
clearly in an expanded view in the lower inset. The upper 
inset shows the variation of fractional discontinuity in TL at 

threshold A ^ t t ] h f = J an( l °f with temperature. 

The two-well design of Fig. [IJ which features an 
intrawell-phonon (IP) depopulation scheme as opposed to 
the RP depopulation scheme of Fig. [5J has also been ex- 



perimentally realized recently [23[ . The design-bias band 
diagram of the realized structure is shown in Fig. [H^a) 
including higher energy parasitic levels that are believed 
to contribute to electron transport at elevated tempera- 
tures. The experimental I-Vs and IZ-Vs from a represen- 
tative ridge laser in cw operation are shown in Fig.[6]^b). 
The temperature variation of the I-Vs shows that J max 
for this IP QCL decreases slightly with temperature up 
to T ~ 45 K and subsequently increases steeply with 
temperature until and beyond the T max cw ~ 63 K. This 
unique dependence of J max with T, also shown explic- 
itly in the inset of Fig. [Hb), is different from that of the 
RP QCLs in Fig. [4] and is attributed to the onset of a 
temperature degradation mechanism that is in addition 
to the reduction in T32 due to thermally activated e-LO 
phonon scattering from 3^2. It is postulated that the 
lifetime of the upper state 3 decreases additionally due 
to absorption of non-equilibrium (hot) LO phonons via 
3' — * (4', 5) scattering [23|]. Any such scattered electrons 
quickly relax back eventually into the injector level of 
a neighboring module due to the short lifetime of the 
parasitic levels 4 and 5 (7-4/ ~ t 5 ~ 0.3 ps). Such a leak- 
age mechanism can be considered approximately as an 
effective reduction in T31 in order to apply the results of 
the 3-level model to explain the behavior in Fig. \&(b), 
which then becomes consistent with calculated I-V be- 
havior in Figs. [21b). Note that the upper laser state in 
the RP designs is spatially isolated from the wide injec- 
tor well(s) and hence from the higher energy parasitic 
levels @. Hence, the aforementioned hot-phonon medi- 
ated leakage mechanism is weaker for the RPTQCLs as 
compared to the IP QCL structure of Fig. [6^a). 

B. Discontinuity in differential resistance at 
threshold 

The onset of lasing in a QCL is characterized by a 
slope-discontinuity in its I-V characteristics [12|. This 
is a manifestation of the increased rate of carrier flow 
due to the stimulated emission process from the upper 
to the lower state, such that the population inversion 
is kept constant beyond threshold. The relative change 
in the differential resistance TL = dV/dl of the device 
at threshold is a useful parameter since it is easily and 
accurately measurable, and is shown below to be directly 
proportional to the value of An t h- We can relate the 
externally applied bias voltage across all the periods of 
the QCL to the energy detuning HAys as 



where, N p is the number of repeated periods in the QCL 
structure, and fv is the fraction of per module voltage 
bias V/ N p that appears as the energy difference between 
levels 1' and 3. fv is typically a slowly varying non-linear 
function of the applied voltage V and is a characteristic 
of a particular design as shown in Fig. [7J For the main 
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result presented in this section, the particular form of fy 
is inconsequential and it is assumed to be a constant in 
the bias range right below and above the threshold. Using 
equation (fT2"|) the differential resistance 1Z becomes 



dV _ HN p dE 1 
dl ~ fv\e\ 



dl 



(13) 
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FIG. 7: Detuning energy as a function of the voltage 

per QCL module (= V/N p in equation 1 12|) for the tightbind- 
ing wavefunctions of the 4-level resonant-phonon design of 
Fig. HJa) and the 3-level intrawell-phonon design of Fig. |6ja) 
respectively. Small circles represent calculated values, which 
are overlaid with straight line fits to the data points close to 
the 1' — 3 resonance. The fits indicate that is linearly 

related to V/N p close to the 1' — 3 resonance (i. e. fv in equa- 
tion[l2]is a constant). Typically, lasing threshold occurs after 
the 1' — 2 resonance (~ 30 mV/module) beyond which the 
current continues to increase up to the 1' — 3 resonance when 

Jmax OCCUrS. 

In equation (|13p. we have dropped the term corre- 
sponding to ^j- since fy is approximately a constant 
as can be seen from Fig. [7] The calculations in Fig. [7| 
are done without accounting for band-bending due to 
dopant-carrier segregation, which is typically insignifi- 
cant due to the low-doping in such QCL structures [l9l |. 
The value of ^j- is likely to be negligible in comparison 
to the leading term of equation (fl~3|) even if band-bending 
were found to be significant for a particular design. We 
can now derive analytical expressions for the ratio as well 
as the fractional change in differential resistance just be- 
low and above (TZ^) threshold using equation ([5|) 



to calculate 



dA, 



dl 



The individual expressions for TZ th 

"-th 



and 



and lZ^ h are complicated; however, the ratio 

hence the fractional change in 1Z at threshold A Jj Lth is 
expressed in a surprisingly concise form as [221 ] 



A^th 




Note that equation (|15|) is derived for the case of unity 
injection efficiency, i. e. Sly 2 ~ in the 3x3 DM equa- 
tion ([2]). Hence, within the approximation that the cur- 
rent flow due to the 1' — > 2 channel is negligible in com- 
parison to the 1' — ► 3 channel at threshold, the expression 
derived for / ^" th holds true irrespective of the nature of 
the RT transport (coherent or incoherent) across the in- 
jector barrier. At low temperatures, r%\ <C T31, and hence 
2A " th gives an absolute measurement of the pop- 

^tot 



ATC t i 



ulaf ion inversion in the laser as a fraction of the combined 
populations of the injector level 1' and the laser levels 3 
and 2 (n to t = ny + n 2 + m). 

The slope efficiency of a QCL's optical power output 
is written as 



dP 



dl 



|e| (a w + Q; m ) 



V 



(16) 



where cuo is the lasing frequency, a m is the radiative (mir- 
ror) loss in the cavity, a w is the material (waveguide) 
loss, and 77 is the internal quantum efficiency of the QCL 
structure, 77 < 1 due to the fact that n.3, and hence the 
non-radiative component of the current / (= \e\n^/T3), 
continues to increase with / above threshold (even as 
An = ri3 — ri2 remains fixed at An t h), which causes the 
current above threshold (/ — 2th) to be not entirely due 
to radiative transitions. An expression for 77 is derived as 



I21 

T32 



where, 



1+121 

T31 



dn 3 

dl lj ='t1 
dn 3 , 
dl 



= 1 



1 1=1:, 



dn 3 I 
dn 3 I 

di 1/=/;, 



1 T21T31 
|e| (t 2 i + T31) 
1 

1 r r 3 



(17) 



(18) 



Equation (jTTJ) for 77 can also be derived from a rate equa- 
tion model since 773 can be expressed in terms of / re- 
gardless of the nature of 1' — > 3 RT transport. In litera- 
ture [H, Hil , A £' th has similarly been derived from a rate 
equations approach, arguing that n.3 oc A1/3, and hence 
the applied voltage V , in which case becomes same 

as 77 (for unity injection efficiency), and hence an indica- 
tor of the ratio of laser level lifetimes T21/Y32. However, 
the expression for 



AU t 



in equation (|15 



„ m ^ 4 uu,uivn is starkly dif- 
ferent from that of 77 in equation (jTTJ) . While 77 — > 1 for 
t 2 i -> 0, but -> ^f-- Also, 77 strongly depends 

on the upper state to lower state lifetime T32, whereas 
A ^ th is independent of T32. Hence, the result derived 
with the 3-level DM model has important implications 
in the way the slope discontinuity is interpreted for the 
I-V of a QCL at threshold. In essence, A ^' th is an indica- 
tor of the population inversion at threshofd and not the 
ratio of level lifetimes, which are two different aspects of 
laser operation. 

By means of Fig. [5J we now show that the assump- 
tion 713 oc Ai'3 which is the basis of the previously held 
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FIG. 8: Illustrative band diagrams for (a) double-barrier 
resonant-tunneling diode structure [30l ] that consists of a sin- 
gle quantum-well sandwiched between degenerately doped 
emitter and collector regions, and (b) a superlattice struc- 
ture [l[ that consists of multiple repeated quantum-wells of 
which only three wells are shown. In (a), the shown Fermi 
distributions (indicated by thin red lines) for the carrier pop- 
ulations apply to a quasi-continuous carrier distribution in 
three-dimensions, whereas in (b), the shown carrier distribu- 
tion applies only to the two-dimensional momentum space in 
the x — y plane (since the carriers are confined in the z direc- 
tion), and is typically Maxwell-Boltzmann like (i. e. the Fermi 
energy [i n lies much below the energy of the bottom of the 
subband E n indicated by thick horizontal lines). In (a), V 
is the voltage applied between the emitter and the collector 
that appears as the difference between the respective Fermi 
energies, and in (b) V represents the voltage per repeated 
module of the superlattice that appears as the difference in 
the energy of the bottom of the subbands. 



belief of A j^ th being equivalent to rj, is incorrect. For 
the conventional resonant-tunneling diode (RTD) struc- 
tures, carriers tunnel from the three-dimensional states 
in the emitter into the quasi-two-dimensional states in 
the well. Considering only first-order tunneling processes 
(i. e. conservation of in-plane momentum in tunneling 
through the barriers), the current flow is proportional to 
the Fermi energy /xe in the emitter and hence the applied 
voltage V (as shown in Fig. Oa), see Ref. HH for exam- 
ple). In such a case, the population n\ of the subband 1 in 
Fig. [3a) will indeed be proportional to V. However, for a 
semiconductor superlattice structure of which QCLs are 
specific examples, the current flow will be negligible if the 
subbands are significantly off-resonance in the neighbor- 
ing quantum-wells even if the carrier distributions over- 
lap in energy. This is because total energy cannot be con- 
served in such a tunneling process since the in-plane mo- 
mentum conservation needs to be satisfied. The applied 
bias field is in the z direction and hence affects the bot- 
tom energy of the subbands and their z-wavef unctions, 
however, the quasi-Fermi levels within the respective sub- 
bands represent the in-plane carrier distribution and are 
not directly affected by the applied voltage V. As a re- 
sult, the subband populations cannot be assumed to be 
directly proportional to V. Instead, the subband pop- 
ulations and their thermal distributions (and hence the 



respective quasi-Fermi levels) are determined by the var- 
ious intersubband and intrasubband scattering processes 
and also by the resonant-tunneling current flow induced 
by the coupling of the subband wavefunctions across the 
barriers as determined by the anticrossing energies 
and detuning energies HAij . Note that higher order cor- 
rections to the resonant-tunneling current flow that relax 
the requirement of in-plane momentum conservation are 
not included in this simplified picture [IH |32| . 
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FIG. 9: Variation of the fractional change in differential re- 
sistance at threshold ^^ tiL with T31 calculated using the ex- 
pression in equation (|15[l with same typical values of different 
parameters as in Fig. [3] 



Figure [9] shows typical variation of / ^* h with T31 ac- 



cording to equation (|15[) . which shows that 



ATCt; 



de- 



creases rapidly as T31 becomes smaller and approaches 
T21. The trend in Fig. [9]could be compared to the exper- 
imental variation of / ^* h with temperature for the two- 
well QCL shown in the upper inset of Fig. H{b), which 
indicates a steady degradation of its effective T31 with 
temperature due to a predicted hot-phonon effect [23l |. 
Also note that </ max and A ^ h in Fig. [6^b) have a corre- 
lated behavior where J max T as !• This validates 
the expression in equation |5]) for operation above thresh- 
old, which suggests J max oc 1 1 



AR th 



where J max is 



the current flowing at 1' — 3 resonance (A1/3 = 0) and 
equation (fTS)) is used for substitution. For this IP QCL 
also, gain bandwidth broadening with temperature is less 
likely to be the dominant temperature degradation mech- 
anism. Otherwise, the value of A J?~ til should increase ac- 
cording to equation (|I5p with the requirement of a larger 
population inversion to reach threshold. 



III. OPTICAL GAIN SPECTRUM OF 
RESONANT-PHONON TERAHERTZ QCLS 

Terahertz QCL designs with RP depopulation typi- 
cally have broad gain lincwidths. The spontaneous emis- 
sion spectra from one of the earliest such designs had 
a - 2 THz full-width half-maximum (FWHM) linewidth 
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4— level resonant-phonon THz QCL design 




FIG. 10: The 4-level resonant-phonon terahertz QCL struc- 
ture reproduced from Fig. |4ja) with a different level num- 
bering scheme for simplicity of presentation in the following 
discussions. The radiative transition is from 4 — ► 3, and de- 
population of the lower level is via 3^2 RT and 2 — > 1 e-LO 
phonon scattering where E21 ~ Swlo- 



for a ~ 4 THz QCL active region [Ti|. Experimental 
spontaneous emission data for some of the more recent 
high performance designs [l(J [HI is not available be- 
cause of the difficulty in observing sub-threshold opti- 
cal signal in metal- metal waveguides [34j], both because 
of the low loss in these waveguides that causes lasing 
soon after the upper level is populated with electrons, 
and also due to the poor out-coupling efficiency of such 
cavities [3S| . However, simultaneous lasing of modes sep- 



arated by frequencies greater than 0.5 THz [20j, l36| sug- 
gests that large gain exists over a broad bandwidth that 
is a significant fraction of the center frequency in a homo- 
geneously grown RPTQCL design (in contrast to hetero- 
geneous cascade designs that have been used to realize 
broad gain bandwidths at mid-infrared frequencies [3?i]). 

In this section, we use density-matrices to numerically 
estimate the optical gain spectrum in a RPTQCL de- 
sign and find reasons for broadening of the gain spec- 
trum. Similar to the I-V calculations in section HH the 
simplified DM calculations are shown to be an effective 
method to understand the role of various design parame- 
ters toward the optical response of a terahertz QCL gain 
medium, as opposed to the nonequilibrium Green's func- 
tion technique that incorporates a full quantum theory 
of gain [H, [3!|. Similar calculations have been used 
recently to determine the optimal injector coupling in 
mid-infrared QCLs [4(|, and also to describe the electri- 
cal transport characteristics of RT extraction based ter- 
ahertz QCL structure [3]. We consider the simplest RP 
design with 4-levels per QCL period as shown in Fig. |TD] 
To compute the gain (or loss) spectrum, the linear re- 
sponse to a sinusoidal electric field perturbation of the 



form E = z£(e 



e ) is calculated by incorpo- 



rating the electric-dipole interaction term for coherent 
coupling of the radiative levels. Only the electric-field 
component in the growth direction z is considered due to 
the intersubband polarization selection rule. The time 
evolution of the DM for the 4-level design is then written 
as 22] 



dt 
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(19) 



Equation (|19[) is written similarly as equation ([2]) with 
additional modifications due to the sinusoidal electric- 
field perturbation. We consider linear response within 
the rotating- wave approximation (i. e. \lj — lo^\ lo) in 
which case the affected coherences can be explicitly writ- 
ten to have a sinusoidal time variation as shown with 
the slowly- varying amplitudes p m n- The radiative lev- 
els 3 and 4 are now coherently coupled through the off- 
diagonal electric-dipole interaction term e|r • E that has 
an amplitude |e|z43^, where Z43 = (4|z|3) is the dipolc- 
matrix element for the said tightbinding levels. Since 
(l',4) and (3,2) are coherently coupled due to non-zero 
f2i'4 and ^32 terms respectively, the coherences corre- 



sponding to (l',3), (l',2), and (4,2) acquire a time- 
harmonic character due to the time-harmonic (4, 3) co- 
herent coupling. The sinusoidal component of the co- 
herences is written explicitly in the ansatz /9(i',4 3 2)5 hi 
which case dp ^ n = in the steady-state for the slowly 
varying amplitudes p mn . Note that P(\\a,3,2) 011 the right- 
side of equation p^|) is the same as the left-side matrix 
appearing with the time derivative. 

We made some simplifying approximations in writing 
equation (fl9|) . Similar to before, for simplicity of calcu- 
lations, we assume a unity injection efficiency (f^i'3 ~ 0), 
also a unity collection efficiency (^42 ~ 0), and neglect 
the backscattering terms rx2 and T34. For non-zero fii'3 
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and/or ^42, equation (|19[) will become more complicated 
since some of the non-harmonic terms in the p^, 432) 
matrix will additionally acquire time-harmonic character 
and similarly some of the time-harmonic coherences will 
additionally acquire a constant value in the steady-state, 
thereby effectively increasing the number of independent 
variables to be solved for the system (currently 16 for 
the written equation). Also, a parameter T41 is included 
to incorporate the effect of any indirect parasitic scatter- 
ing channels from 4 to 1, since direct 4 — ► 1 transport 
is otherwise only allowed via inter-sub-modulc tunneling 
across the collector barrier for the chosen set of tight- 
binding basis functions. We can now solve this equation, 
which, in addition to yielding values of current and level 
populations, also yields the electrical polarization iV(t) 
induced due to the external optical field E as 

V = e oX (uj)£e+^ t + c.c. 
_ ~\e\(z) 



V ac 
je|f43 



{[P43(w) + P34(-w)] < 



•} (20) 



where V ac is the volume of the active region. The induced 
electrical susceptibility x(w) becomes 



-|e|z43 



V ac e £ 



043 M +p34.(-u)} 



(21) 



which is independent of the amplitude £ for small values. 
The optical gain coefficient g(u>) (in meter -1 ) is related 
to the imaginary part of the susceptibility x"( w ) as 



'.x")<"? 



n r c 



(22) 



where n t is the refractive index of the medium. In the 
following, we show results of the computed gain spectrum 
as a function of several different parameters of interest. 



A. Optical gain spectrum as a function of injector 
anticrossing 

Figure [TT] shows the numerically computed gain spec- 
tra for the 4-level RPTQCL design as modeled by 
equation (fT9|) for different values of injector anticross- 
ing 2Ml\n. For the 4 — > 3 radiative transition, the 
FWHM frequency linewidth due to scattering is given 



by A^ s , 



1 

2T4 



1 



which yields a value 



_ 2^3, off 

of 0.67 THz (/iAi/ scatt ~ 2.75 meV) for the chosen pa- 
rameters. The value of T 2 * = 0.75 ps was assumed for 
the above value of the linewidth to approximately agree 
with those observed typically for bound-to-continuum 
THz QCLs [HI, HH, HH whose linewidths are more likely 
to be scattering limited. As 2K£lvi becomes as large 
as /iA^ sca tt, injector transport becomes more coherent 
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FIG. 11: Computed gain spectra for the 4-level structure of 
Fig.[l0]for different injector anticrossings 2Kl 1 i^, evaluated at 
1' — 4 resonance (A1/4 = 0). Levels 3 and 2 are also taken to be 
at resonance (A32 = 0) and a small value of 2MI32 = 2 meV is 
chosen to limit additional broadening due to a coherent 3^2 
RT process (as shown in Fig. ll2b h All other parameters have 
same equivalent values as in Fig. except T21 = 0.2 ps that 

results in an effective lower state lifetime Tz, e « {= ^f^ij ~ 
1.0 ps, which can also be estimated from the equivalent of 
equation [fTT)l. Also, £43 = 16.5 meV (4 THz). The values of 
FWHM linewidth Az/, and population inversion An43 = 714 — 
723 obtained as a result of the calculation are also indicated 
alongside each of the curves. 



and the gain spectrum becomes additionally broadened 
due to anticrossing splitting of the linewidth as can be 
seen from Fig. [TTJ Also note that the population inver- 
sion An approaches the maximum value given by equa- 
tion (J7]) for coherent injection, which, for the present case 



yields An„ 



/ 1-T 3 , c ff/T4 3 \ 

V2+r3. ctf /T4 3 y ' hot 



0.28 n tot . The 



v 2+T3 i6 ff /Y43 f 

area under the gain curve is proportional to the popula- 
tion inversion Ari43, therefore, any additional broaden- 
ing decreases the peak gain. This suggests that a value 
of 2M7]/4 much larger than hAv sca tt might in fact di- 
minish the performance of a design. However, as will be 
shown in the next section, RPTQCL designs anyway have 
broad linewidths due to a coherent RT assisted depopu- 
lation process, in which case having 2hflin > hAv sca t t 
is less likely to cause additional broadening of the gain 
linewidth. Consequently, the injector anticrossing value 
for the RPTQCL designs is chosen based on other design 
parameters such as the 1' — > 3 parasitic current leak- 
age [Io| rather than the concerns about gain broadening. 



B. Optical gain spectrum as a function of collector 
anticrossing 

Figures [T^f a) and[T2Tb) show the computed gain spec- 
tra for different values of collector anticrossings 2hfls2, 
and energy detuning fiA 32 , respectively. The best RP- 
TQCL designs typically have collector anticrossing values 
in the range of 4 — 5 meV. This is to maintain a short 
effective lower state lifetime so as to maximize J max / Jth, 
since J max oc T[~^ er for coherent injection (equation[5]), al- 
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FIG. 12: Computed gain spectra for the 4-level QCL design 
of Fig. \W\ for (a) different values of collector anticrossings 
2HQ32 at resonance (A32 = 0), and for (b) different values 
of detuning hAs2 at 2HQ32 = 5 meV. Levels 1' and 4 are 
assumed to be at resonance (A]/ 4 = 0) and a small value of 
2hQ 1 i 4 — 1.5 meV is chosen to limit additional broadening 
due to a coherent 1' — * 4 RT process (as shown in Fig. IX 1 p i . 
Other parameters are the same as in Fig. 1111 The values 
of FWHM linewidth Av, effective lifetime of the lower state 

T3 iC fr nf T21 ) ' ano - P P u l a ti on inversion A7143 = 714 — 713 ob- 
tained as a result of the calculation are also indicated along- 
side each of the curves. 



though at the cost of a broader linewidth. For example, 
the 2hn 32 = 5 meV spectrum in Fig.H^a) has a FWHM 
linewidth of 1.8 THz, a value almost twice as large as 
the scattering linewidth Au sca tt ~ 1 THz (calculated us- 
ing T3 je ff ~ 0.3 ps). This shows that the anticrossing 
splitting of the gain spectrum due to the coherent RT 
assisted depopulation process is the main cause of the 
broad linewidths associated with RPTQCLs [H, H3, H| • 
As seen from Fig. [T27 b). note that linewidth broadening 
happens only close to the 3 — 2 resonance at I1A32 ~ 0. 
For |fiA 32 | > and |ftA 32 | ^ £43 , the width of the spec- 
trum remains narrow, since the 4^2 radiative transition 
(which is indirect, and due to 4 <-> 3 and 3^2 coher- 
ent coupling) is considerably detuned from the stronger 
4^3 radiative transition (which is direct, and due to 
4^3 coherent coupling), and hence contributes neg- 
ligibly to the gain. For a coherent 3 — > 2 RT process 
(i. e. large SI32), the peak gain occurs at a frequency 



Wp C ak > E i3 /h for HA 32 < and w pcak < E i3 /h for 
HA 3 2 > 0. The frequency shift of w pca k away from E± 3 /h 
is due to the energy splitting due to the 3 — 2 anticrossing 
close to resonance. 
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FIG. 13: (a) Computed I-V , and gain spectra at selected bias 
points for the 4-level QCL in Ref. [nj. These results are ob- 
tained by solving equation (|19|l in the steady-state (-& — ► 0), 

whereby, I = |e| ^fj^ + fjj) an d the gain is calculated using 

equation (|22|) . A value of g t h = 40 cm" 1 is assumed similar 
to that for Fig. [3] [2^|, following which the stimulated life- 
time 7~st is determined to keep the peak gain constant as gth 
beyond threshold. Non-radiative lifetimes are calculated by 
including e-LO phonon scattering only whereby an electron 
temperature of 100 K is assumed in the upper level 4. Also, 
T 2 * ~ 0.75 ps, and fi^ « Q42 ~ 0. (b) Experimental I-V, 
and spectra (shown as insets) at selected bias points measured 
at a heat-sink temperature of 20 K in cw mode. The device is 
a 100 /im x 1.39 mm metal-metal waveguide ridge laser from 
the same active region as in Ref. fxcL The higher than ex- 
pected bias voltages for the experimental I-V are likely due 
to additional voltage drop at the electrical contacts. 

We now show comput ed g ain spectra of the 4-level RP- 
TQCL design from Ref. lfj in Fig. fl~37 a) as a function of 
applied electrical bias. For a comparison, experimentally 
measured cw I-V and lasing spectra of a ridge laser pro- 
cessed from the same active medium is also shown in 
Fig- fTBT b). For the calculation in Fig. [TBT a). the radia- 
tive energy £43 varies linearly from 12.4 meV (3.0 THz) 
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at 8 V to 15.7 meV (3.8 THz) at the 1'- 4 resonance bias 
of 12 V due to the large value of the Stark shift in a diag- 
onal radiative transition. However, the peak gain occurs 
close to a frequency of 4 THz at almost all bias condi- 
tions for the calculated spectra. This behavior is due to 
the coherent nature of the depopulation process. The 
3 — 2 energy detuning /1A32 also varies linearly with the 
applied bias, from -4.7 meV at 8 V to 2.4 meV at 12 V. 
Consequently, as shown by the calculations in Fig. [T27 b). 
the gain peak is pushed to higher frequencies than the ra- 
diative separation S43 at low-bias that effectively masks 
the Stark shift in E43. This is affirmed from the las- 
ing spectra in Fig. [T37 b) that shows that the ^3.9 THz 
mode is excited at all bias. In contrast, the Stark shift 
is clearly emphasized in the lasing mode spectrum of the 
diagonal bound-to-continuum terahertz QCLs [7j due to 
the different depopulation mechanism in such designs. 

As an evidence of the broad gain linewidth in this QCL, 
at some middle-bias points, additional lasing spectra are 
excited close to a frequency of 3.3 THz that leads to si- 
multaneous cw lasing of modes separated by ~ 0.6 THz. 
This behavior could be attributed to the double-peaked 
gain spectra that emerges close to the 3 — 2 resonance 
(A 32 ~ 0), which, for this design, happens at the middle- 
bias points. It may be noted that we have observed sim- 
ilar dual-frequency lasing behavior in some versions of 
the 5-level RPTQCL design of Fig.^b). Whether or not 
this happens depends on the relative alignment of the 
injector and collector anticrossings as a function of the 
applied electrical bias. 

We should note that recently a much smaller FWHM 
linewidth of ~ 0.6 THz has been measured for a similar 
4-level RPTQCL design Q using a THz time-domain 
spectroscopy technique, which could be due to the follow- 
ing reasons. The QCL structure in Ref. has a smaller 
collector anticrossing 2hQ 32 ^3.7 meV as opposed to a 
value of 2^32 ~ 4.7 meV for the design in Ref. [lOj that is 
analyzed in Fig. Q1J] Also, it is possible for a given QCL 
structure that the bias range of lasing operation may not 
sweep through its collector resonance (A32 = 0) depend- 
ing on the relative alignment of the injector and the col- 
lector anticrossings for the grown structure, in which case 
the gain spectra may not show additional broadening due 



to the 3 — > 2 RT as seen from some specific calculations 
for non-zero SA32 in Fig. [T27b) . We would also like to 
note that the present calculations were done within a 
rotating-wave approximation (\lu — uj^I <C uj), which be- 
comes less accurate once the gain linewidth becomes a 
significant fraction of the center frequency. Hence, the 
true linewidths may be somewhat narrower than those 
calculated in this section. Nevertheless, within the as- 
sumptions considered, the close tracking of the calculated 
and experimental spectral characteristics of the 4-level 
RPTQCL design in Fig. Q]|] establishes the importance of 
the relatively simple density-matrix model developed in 
this section for estimating the optical gain spectrum of a 
terahertz QCL structure. 

In conclusion, we have presented simplified density- 
matrix transport models to describe resonant-tunneling 
transport in terahertz QCLs. Due to the closely spaced 
energy levels in terahertz QCL structures, coherence 
plays an important role in the resonant-tunneling mech- 
anism, which is incorporated well within the presented 
framework. A 3-level model was developed to derive cur- 
rent transport through the injector barrier of any gen- 
eral QCL design (which applies to mid-infrared QCLs as 
well). Useful expressions were derived for current flow 
above and below threshold that could directly be used to 
analyze the experimental behavior of some representa- 
tive QCL designs. Based on experimental observations, 
we have been able to speculate about some of the dom- 
inant temperature degradation mechanisms in phonon- 
depopulated terahertz QCL designs. We have extended 
the density-matrix model to estimate the gain spectra of 
resonant-phonon terahertz QCLs. A coherent resonant- 
tunneling assisted depopulation process is identified to be 
the primary cause of the broad gain bandwidths typically 
observed in such QCLs. 
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